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ABSTRACT 

In this paper we report 23 magnetic field measurements of the B3IV star HD 23478: 
12 obtained from high resolution Stokes V spectra using the ESPaDOnS (CFHT) 
and Narval (TBL) spectropolarimeters, and 11 from medium resolution Stokes V 
spectra obtained with the DimaPol spectropolarimeter (DAO). HD23478 was one of 
two rapidly rotating stars identified as potential “centrifugal magnetosphere” hosts 
based on IR observations from the Apache Point Observatory Galactic Evolution Ex¬ 
periment survey. We derive basic physical properties of this star including its mass 
(M = Mq), effective temperature (Te® = 20±2kK), radius {R = 2.7^1'^ Rq), 

and age (rage = Myr). We repeatedly detect weakly-variable Zeeman signatures in 
metal, He and H lines in all our observations corresponding to a longitudinal magnetic 
field of (Bz) « —2.0kG. The rotational period is inferred from Hipparcos photome¬ 
try (Prot = 1.0498(4) d). Under the assumption of the Oblique Rotator Model, our 
obsevations yield a surface dipole magnetic field of strength By ^ 9.5 kG that is ap¬ 
proximately aligned with the stellar rotation axis. We confirm the presence of strong 
and broad Ha emission and gauge the volume of this star’s centrifugal magnetosphere 
to be consistent with those of other Ha emitting centrifugal magnetosphere stars based 
on the large inferred Alfven to Kepler radius ratio. 

Key words: Stars: early-type. Stars: magnetic fields, Stars: individual: HD 23478 


* Based on spectropolarimetric observations obtained at the 
Canada-France-Hawaii Telescope (CFHT) which is operated by 
the National Research Council of Canada, the Institut National 
des Sciences de I’Univers (INSU) of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii, 
observations obtained using the Narval spectropolarimeter at the 
Observatoire du Pic du Midi (France), which is operated by the 


1 INTRODUCTION 

Upon the first detection of its magnetic field, the strong, 
broad, and variable Ha emission line of the main sequence 

INSU, and observations obtained at the Dominion Astrophysical 
Observatory, NRC Herzberg, Programs in Astronomy and Astro¬ 
physics, National Research Council of Canada. 
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Table 1. Observations of HD 23478 obtained with ESPaDOnS and Narval. The phases are calculated using Eqn. Q. SNRs per 1.8 km s“^ 
spectral pixel are reported at 5400 A. The three rightmost columns list the longitudinal magnetic field measurements obtained from H/3 
and LSD profiles generated from a He+metal line mask along with the detection status according to the criteria of |Donati et al.| l |1997[ |: 
definite detection (DD), marginal detection (MD), and no detection (ND) (see Section]^ 


HJD 

Phase 

Total Exp. 
Time [s] 

SNR 

[pix-i] 

Instrument 

(-Bz>H/3 

[kO] 

{Bz ) He+metal 

[kG] 

Detection 

Status 

2456884.116 

0.304 

560 

473 

ESPaDOnS 

-1.52 ± 0.27 

-1.98 ±0.25 

DD 

2456908.146 

0.194 

560 

500 

ESPaDOnS 

-1.68 ±0.24 

-2.15 ±0.22 

DD 

2456909.134 

0.135 

560 

381 

ESPaDOnS 

-1.42 ± 0.32 

-2.33 ±0.29 

DD 

2456909.144 

0.145 

560 

346 

ESPaDOnS 

-1.89 ±0.29 

-2.26 ±0.26 

DD 

2456925.512 

0.736 

1200 

402 

Narval 

-2.19 ±0.37 

-1.98 ±0.36 

MD 

2456925.691 

0.907 

1200 

291 

Narval 

-2.32 ± 0.28 

-2.49 ±0.30 

DD 

2456962.025 

0.517 

560 

409 

ESPaDOnS 

-0.95 ± 0.34 

-1.93 ±0.26 

DD 

2456968.937 

0.101 

560 

409 

ESPaDOnS 

-2.13 ±0.31 

-2.44 ± 0.27 

DD 

2456972.833 

0.813 

560 

353 

ESPaDOnS 

-1.94 ±0.39 

-2.52 ±0.37 

MD 

2456972.843 

0.822 

560 

352 

ESPaDOnS 

-2.98 ±0.41 

-2.84 ± 0.37 

MD 

2456973.099 

0.066 

560 

497 

ESPaDOnS 

-2.05 ± 0.27 

-2.11 ± 0.22 

DD 

2457034.695 

0.740 

560 

461 

ESPaDOnS 

-1.63 ±0.29 

-1.99 ±0.27 

MD 


B2 star cr Ori E was suggested to be a natural consequence 
of a plasma co-rotating with the star well beyond its surface 
( Landstreet &: Borra|1978 l. Since then, a significant number 
(> 13) of magnetic mid- to early-B stars (i.e. Tea < 25 kK) 


exhibiting Ha emission were discovered (Brown, Shore & 
Sonneborn|| 1985 Shore et al.||199(5 |. These stars are of par¬ 
ticular interest because they provide a unique means of prob¬ 
ing the mass-loss and stellar winds of B-type stars. While 
the strength of the magnetic field likely plays an important 
role in determining if a magnetic B star exhibits Ha emis¬ 
sion, the key distinction between the stars with emission 
and those without is rotation - the emission line stars are 


all rapid rotators having rotational periods < 1.5d (Petit 


et al. 12013 1. 


These correlations between rotation and the presence of 
Ha emission can be understood in a general framework in 
which the magnetospheres of massive stars are classified as 
either dynamic magnetospheres (DM) or centrifugal magne¬ 
tospheres (CM) based on the ability of the field to confine 
the wind and the degree of criticality of the stellar rotation 


(ud-Doula & Owocki 2002 Petit et al. I 20131. Petit et al. 


(20131, and more recently Shultz et al. ( |2014 (, demonstrate 
that stars like a Ori E are those with the largest magne- 
tospheric volumes: those with the largest Alfven radii com¬ 


bined with the smallest Kepler radii. Shultz et al. (20141 


estimated that Ha emission is never observed for typical 
magnetic B stars with Ra/Rk < 10, while it is frequently 
observed in the spectra of stars with ratios above this thresh¬ 
old. 


Recently, Eikenberry et al. (20141 identified two emis¬ 


sion line B stars from near-infrared (nIR) spectra obtained 
in the context of the Apache Point Observatory Galactic 
Evolution Experiment (APOGEE) survey: HD 345439 and 
HD 23478. They found that both stars exhibited strong and 
broad emission at all nIR H lines along with relatively large 
projected rotational velocities of « 270 and 125 km s“^, re¬ 
spectively. Although lacking any information about their 
magnetic fields, Eikenberry et al. (20141 concluded that the 


sonable evidence for a magnetospheric origin of the observed 
nIR emission. 

The principal aim of the present study is to search for 
the presence of a magnetic field in the brightest of these two 


similarities between these stars and a Ori E provided rea- 


stars (HD 23478), and to test the hypothesis of Eikenberry 
|et al.| ( |2014[ ) that it is similar to cr Ori E in that it also be¬ 
longs to the more general class of Ha emitting GM stars. In 
Section we describe the spectropolarimetric observations 
used in this study. In Section we estimate the physical 
parameters of HD 23478 using photometric data from the 
literature along with newly-acquired high resolution spec¬ 
troscopy. We then analyze in Section variability in Hip- 
parcos epoch photometry in an attempt to infer the star’s 
rotational period. In Section we compute the longitudinal 
magnetic field and search for its variability. Section dis¬ 
cusses the spectral line emission and variability. In Section 
[^we derive the stellar magnetospheric parameters. Einally, 
in Section we evaluate our hypothesis for a CM origin of 
HD 23478’s emission. 


2 OBSERVATIONS 

2.1 High resolution spectropolarimetry 

HD 23478 was observed during eight nights between Aug. 
2014 and Jan. 2015 using the ESPaDOnS and Narval spec- 
tropolarimeters installed at the Canada-France-Hawaii Tele¬ 
scope (CFHT) and Telescope Bernard Lyot (TBL), respec¬ 
tively. Both instruments have a resolving power of R ~ 
65 000 acquiring spectra in circularly polarized light span¬ 
ning the visible range of 3 600 — 10 000 A. The Heliocentric 
Julian Dates (HJDs), exposure times, and signal-to-noise ra¬ 
tios (SNRs) of these observations are listed in Table 

Each circular polarization observation consists of four 
subexposures which the Libre-ESpRIT pipeline ( [Donati 
et al.||19^ I automatically reduces yielding the final Stokes I 
and V spectra. The majority of the observations (10 out of 
12) were obtained with ESPaDOnS using a total exposure 
time of 560 s. The remaining two observations were obtained 
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Table 2. Observations of HD 23478 obtained with dimaPol. The 
phases are calculated using Eqn. 0 and the magnetic field mea¬ 
surements correspond to the longitudinal field derived from H/3 
Stokes V observations. 


HJD 

Phase 

Total Exp. 
Time [s] 

SNR 

[pix-i] 

{-Bz}h/3 

[kC] 

2456680.700 

0.537 

5400 

300 

-1.78 ±0.56 

2456681.630 

0.422 

3000 

250 

-0.92 ±0.72 

2456963.734 

0.146 

3000 

440 

-1.62 ±0.27 

2456971.888 

0.911 

4800 

400 

-2.88 ±0.36 

2456972.802 

0.782 

4800 

520 

-1.82 ±0.26 

2456973.758 

0.693 

4800 

430 

-1.29 ±0.41 

2456974.846 

0.731 

4800 

410 

-1.37 ± 0.38 

2456976.032 

0.859 

4800 

400 

-2.08 ±0.29 

2456991.962 

0.033 

4800 

480 

-1.95 ±0.30 

2456992.711 

0.748 

4800 

310 

-1.86 ±0.49 

2457085.723 

0.347 

4800 

550 

-1.11 ± 0.33 


with Narval during a single night using a total exposure time 
of 1200 s. 


2.2 Medium resolution H/3 spectropolarimetry 


Eleven spectropolarimetric observations were also obtained 
using the dimaPol medium resolution spectropolarimeter lo¬ 
cated at the Dominion Astronomical Observatory (DAO) 
from Jan. 2014 to March 2015. The instrument has a re¬ 
solving power of 7? ~ 10, 000 and is optimized to acquire 


spectra over a wavelength range of 4 700 — 5 300A (Monin 
let al.|2012[ ). The HJDs, exposure times, SNRs, and longitu¬ 
dinal magnetic field measurements derived from H/3 Stokes 
V profiles are listed in Table 0 


3 PHYSICAL PARAMETERS 
3.1 SED fitting 

Photometry spanning a wide range of wavelengths is avail¬ 
able for HD 23478 thereby providing a means of constrain¬ 
ing both the effective temperature and surface gravity. These 


observations include monochromatic Johnson UBV | Ander- 
|son fc Francis|2012[), 2MASS JHKs (|Cutri et al.|2003[ ), and 


TDl ultraviolet ( |Thompson et al.||1978 l photometric mea¬ 

surements. Short (SWP imager: 1150 — 1975 A) and long 
(LWP imager: 1900 — 3125 A) wavelength ultraviolet broad¬ 
band measurements were also recorded by the International 
Ultraviolet Explorer (lUE). The high resolution {5\ = 0.1 A) 
lUE data were rebinned at 5A = 3 A using the prescrip¬ 


tion of Solano (19981 in order to increase the SNR per 


pixel. All of the observed magnitudes were converted to 
physical fluxes using published zero points for Johnson and 


2MASS filters (Bessell 1979 Cohen, Wheaton & Megeath 


2003); both TDl and lUE measurements were obtained 
in the appropriate units of ergss“^cm“^A . Deredden¬ 
ing was applied using the method of |Cardelli, Clayton fc| 
Mathis (19891 using a typical total-to-selective extinction of 


R{V) = A{V)/E{B -V) = 3.1. 

The photometric observations were compared with the 
non-local thermodynamic equilibrium (NLTE) TLUSTY 



3.2 3.4 3.6 3.8 4.0 4.2 4.4 

log(A) [A] 



Figure 1. Top: Comparing the lUE spectra (solid blue) and 
monochromatic fluxes (red points) with the best-fitting TLUSTY 
model (dashed black) convolved by a Gaussian profile to match 
the rebinned lUE resolution. Bottom: Theoretical U, B, and V 
fluxes (black crosses) were calculated by convolving the unbroad¬ 
ened TLUSTY models (solid black) with the filter’s associated 
transmission functions. 



1500 2000 2500 3000 

Wavelength [A] 


Figure 2. Comparisons between the lUE LWP spectrum (solid 
blue) and model SEDs (dashed black) at various effective temper¬ 
atures and colour excesses. Both the models and the observations 
are normalized to the V flux (Fy). 


BSTAR2006 synthetic energy distribution (SED) grid ( |Lanz| 
fc Hubeny|2007 1 for effective temperatures Tefi ranging from 


15 — 30kK and surface gravities log// (cgs) from 3 — 4.75. 
The published grid has Tes and log// resolutions of 1000K 
and 0.25, respectively, that were linearly interpolated to a 
finer grid of 250 K and 0.05 resolution. Solar abundances and 
a 2 km s“^ turbulence velocity were assumed for the model 
grid. The NLTE model SEDs have a resolution > 0.002 A 
- much higher than the rebinned lUE observations. We 
reduced the resolution of the model spectra by convolv¬ 
ing with a Gaussian profile assuming a resolving power of 
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Table 3. Stellar parameters of HD 23478. 


Sp. Typei 

B3IV 

TT [mas]^ 

4.99 ±0.62 

d [pc] 

200+29 

Photometric 


V [magj^P 

6.67t9;°9 

E{B - V) [mag] 

U.oU_0.04 

BC [mag] 

— 1 97 + 0-24 
- 0.22 

Mv [mag] 

n 77 + 0-26 
O'-' '- 0.36 

Afbol [mag] 

—2 74+0-51 
^•'^-0.58 

Physical 


Teff [kK] 

20 ± 2 

log(9) [cffs] 

4.2 ±0.2 

V sin i [km s~^] 

140 ± 10 

logL/Lg 

3.0 ±0.2 

M/Mq 

6 1 + 0 -^ 
0 * 1 - 0.7 

R/Rq 

9 7 + 1.6 

^•'- 0.9 

Tage [Myr] 

3lf 

Trot [d] 

1.0498 ±0.0004 

Vr [km s~^] 

17 ± 5 


Table r eferences: ^ Blaauw &: van Albada ([1963|^ van Leeuwen] 
| |2007| , ^Syfert, Hardie fc Grenchik| ( [1960 , ^Harris ( |1956^ . 


R = (1150A)/(3A) « 380 to match the 3 A lUE bin-width 
at A = 1150 A. Theoretical monochromatic fluxes were cal¬ 
culated by convolving Johnson UBV and 2MASS JHKs 


transmission functions (Landolt & Uomoto 2007 Cohen, 


Wheaton & Megeath 20031 with the unbroadened model 


spectra. 

The best fitting NLTE model, shown in Fig. was 
found to have an effective temperature, surface gravity, and 
colour excess of Tes = 22.00kK, logp = 4.60 (cgs), and 
E{B — V) = 0.35 yielding a reduced value of 1.70. The 
monochromatic fluxes were found to be adequately repro¬ 
duced by a wide range of Tefi, loggr, and E{B — V) val¬ 
ues. Changes in the surface gravity of each model SED were 
found to primarily affect the depths of the absorption lines, 
resulting in minimal changes in the calculated monochro¬ 
matic fluxes. The lUE spectra provided the greatest con¬ 
straint on the fitting parameters. In particular, we found 
relatively large discrepancies between the lUE LWP spectra 
and the model SEDs having Teg < 20 kK and Teg > 24 kK 
at 2,000 J; A < 2, 500 A. This is clearly a consequence of the 
respectively lower and higher values of the required extinc¬ 
tion associated with these models. In Fig. we show the 
LWP spectrum dereddened in order to fit models of 18, 20, 
22, and 24 kK. The 2,000 — 2, 500 A region shows growing 
discrepancies for the lower and higher temperatures. This 
leads us to prefer Teg = 22 kK based on the SED fitting. 


3.2 Spectral line fitting 

Spectral line fitting can provide a further constraint on Teg 
and logg through comparisons between the obtained ES- 
PaDOnS and Narval spectra and synthetic model grids. We 
found that HD 23478’s metallic line spectrum was not suit- 




Figure 3. Comparisons between two DIBs found in the spectra 
of HD 23478 and Oph. 


able for tuning the atmospheric parameters due to a num¬ 
ber of characteristics. The strongest metal lines have shallow 
depths of less than 10% of the continuum with only 3 lines 
deeper than 5%. This difficulty was alleviated to some ex¬ 
tent by using averaged spectra thereby increasing the SNR 
by a factor of « 3.5. 

We also noticed that many line profiles were distorted. 
Sharp asymmetric features having depths < 10% of the 
continuum are prevalent in many spectral orders (Fig. |^. 
These were identified as diffuse interstellar bands (DIBs) by 
conferring with Hobbs et al. 12008 20091 and by directly 
comparing the spectra of HD 23478 and C, Oph, a star which 


is well known to exhibit strong DIB features (e.g. Walker, 
Bohlender fc Krelowski|2000 1. 


In addition to the DIB features, many metal lines show 
peculiar strengths and distorted profiles, likely due to non¬ 
solar surface abundances and nonuniform distributions of 
the abundances of these elements. As a consequence, mod¬ 
elling the metallic spectrum to constrain the atmospheric 
parameters would require a detailed analysis which is out¬ 
side the scope of this work. Therefore we opted to model the 
Balmer lines to provide further constraint on Teg and logg 
because they are more weakly influenced by abundance non¬ 
uniformities. Specifically, Hy and HJ were chosen because of 
their relative insensitivity to emission (as compared to Ha) 
and because they are not directly adjacent to other strong 
absorption features thus reducing normalization errors. 

We relied primarily upon spectral models calculated us- 


ing SYNPLOT, the IDL wrapper for SYNSPEC (Hubeny 


& Lanz 


20111, along with NLTE TLUSTY model atmo¬ 


spheres (Lanz & Hubeny 20071. The grid was linearly in¬ 


terpolated such that a finer resolution of ST^s = 250 K 
and Jlogg = 0.1 (cgs) was obtained. The synthetic spec¬ 
tra were then convolved with a Gaussian profile associated 
with the characteristic R = 65,000 resolving power of the 
ESPaDOnS and Narval spectropolarimeters, corresponding 
to a velocity resolution of Sv ~ 4.6 km s“^. It is noted that 
changes in the He abundance can affect the wings of the 
Balmer lines ( Leone fc Manfre|1997 1. When calculating the 
NLTE H line profiles, we use the minimum He abundance 
of [He/H] = —0.58 required to fit a sample of averaged He i 


lines (where [Hb/HJq = —1.01; Grevesse, Noels & Sauval 


19961. 


As shown in Fig. the best effective temperature ac¬ 
cording to the SED fitting (22 kK) reproduces the Balmer 
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Figure 4. Comparisons between TLUSTY models (blue) and the averaged observed spectrum (black) for H(5 (left) and H 7 (right). The 
bottom frames show the residuals between the observations and the synthetic spectra (Fobs ~ -^syn) 


1,00 


1.00 

^ Cll: 4267 A 




\\ v:5lni = 140 

5 0,95 


3 0.98 

\\ vsini = 150 .77 

Ll 


Ll 

\\ // 

X) 


-Q 

\\ / / 

0,90 


N 

\\\. Jj 



'o 0.96 

\ \V^ \l i 

e 


E 

VA n / 





z 0,85 

■ Hel; 4471 A| // 

-z. 

j/ 


“ .'U 1 11 

0.94 

'vo. 


vsini ^ 140 \ 1 



0.80 

. vsinl = 1 50 \J/ 




4460 4465 4470 4475 4480 4485 4265 4266 4267 4268 4269 

Wovelength [A] Wovelength [A] 






Figure 5. = 22 kK and logg = 4.4 TLUSTY models with a range of nsini values compared with the averaged observed spectrum 

(black). The abundances of each element were adjusted in order to yield the best overall fit to the synthetic spectra. DIBs in both the 
Si II A5056 and Fe II A5169 lines are visible at 5055 A and 5170.5 A. 


line wings reasonably well with loggr = 4.4. However, the 
core depths are seriously underestimated by the model. By 
reducing good agreement between the observed and 
computed wings and core is achieved for Teg = 18 kK and 
logg = 4.0. 


We therefore find that the Balmer lines and SED yield 
different solutions to the atmospheric parameters. Based on 
our analysis, we conclude that the effective temperature of 
HD 23478 is likely in the range 18 — 22 kK and its surface 
gravity between 4.0 and 4.5. Further refinement of these 
parameters will require a more detailed future analysis. 

The colour excess of E{B — V) = O.SOIqq® inferred for 
the Teff = 20 ± 2 kK range is in agreement with both the 
0.28 and 0.34 values reported (without consideration of un¬ 


certainties) by Voshchinnikov et al. (20121 and Beeckmans & 


Hubert-Delplace (19801. For T^g = 22 kK however, a slightly 


higher E{B — V) — 0.35 is required. As was expected, Teir 
was indeed found to be closely coupled with the colour ex¬ 
cess. 


3.8.1 Rotational broadening 


To estimate the rotational broadening, we compared ob¬ 
served phase-averaged He and metal profiles with synthetic 
profiles computed assuming Tefi = 22 kK. We adjusted the 
projected rotational velocity, usini, as well as the abun¬ 
dance of each element in order to obtain a reasonable fit 
to the line profiles. Generally, a usinf of 140 ± 10 km s“^ 
yielded the best overall fit to the line profiles (consistent 


with that of Eikenberry et al. 20141. Fits to selected lines 
for several usini values are shown in Fig. The inclusion 
of either macroturbulent or microturbulent broadening did 
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not significantly improve the fits. The distorted shapes of 
many line profiles are apparent. We also note the remark¬ 
ably poor fit to the He l A4471 line. This discrepancy cannot 
be relieved by changing or logg. It is reproduced in a 
number of other He l lines. We suspect that it is a conse¬ 
quence of stratification of He in the atmosphere of this star 
(e.g. Bohlend^|1989 Farthmann et al. 19941. 


3.3 Hertzsprung-Russell Diagram 


The mass and age of HD 23478 were estimated by comparing 
its location on the Hertzsprung Russell diagram (HRD) with 
theoretical isochrones and evolutionary tracks. 

The luminosity of HD 23478 required to determine its 
location on the HRD was calculated using the de-reddened 
visual magnitude of V = 5.62lQ g 7 mag where the color ex¬ 
cess of E{B — V) — 0.30 discussed in Section [3.2| was used. 
The calculation of the absolute visual magnitude, My, re¬ 
quired a bolometric correction which was approximated by 
the temperature and surface gravity relation by |Nieva] ( |2013[ ) 
applicable to stars with 15.8 ^ Teff ^ 34.0 kK. The luminos¬ 
ity was then found to be L = 1 . 01 q '4 ^ 10® Lq where the 
primary contributions to the uncertainty are from the tem¬ 
perature (through the BC calibration) and parallax (i.e. the 
distance). The temperature of Tes = 20 ± 2kK yielded by 
the SED fitting was used in the final HRD placement. 

The comparisons used Geneva model grids of isochrones 
ranging in stellar age from 3.16 Myr to 12.9 Gyr generated 
by [Ekstrom et al.| | |2012[ ). The associated grid of evolu¬ 
tion tracks contains models ranging in mass from 0.8 to 
120 Mq ; both the isochrone and evolution grids assume a 
solar metallicity of Z = 0.014. An identical resolution, 
Z = Zq grid which included rotational effects for a star 


having u/ucrit = 0.4 (Georgy et al. 20131 was also used in 


the analysis. However, the comparisons between the derived 
mass and radius of HD 23478 using both the rotating and 
non-rotating model grids yielded negligible differences of less 
than 5%. 

The position of HD 23478 on the HRD is shown in Fig. 
where the nearby stellar evolution tracks are shown as 
dotted black lines and the main sequence masses are indi¬ 
cated by black ‘plus’ symbols. The isochrones (indicated by 
red dotted lines), logL, and logTefi imply a range in the 
estimated age and radius for the star of r = 31®^ Myr and 
R = 2.81q2,Rq, respectively. Likewise, the main sequence 
mass of HD 23478 is estimated to be M = Mq. A 

similar radius of R — 2. 7^0 9 , Rq is derived from L and Tgfi 
using the Stefan-Boltzmann law; we adopt this value be¬ 
cause of its larger uncertainty. We note that these mass and 
radius estimates imply a surface gravity of logg = 4.31 q' 4 
which is consistent with the value derived from the spectral 
line fitting procedure described in Section [3.2[ 


4 ROTATIONAL PERIOD 


The Hipparcos Epoch Photometry catalogue (Perryman 
et al.|1997 1 contains 75 measurements of HD 23478 spanning 


a time period of three years with minimum and maximum 
time separations of approximately 20 minutes and 174 days, 
respectively. The measurements are therefore sufficient for 
detecting photometric variability across a wide range of time 



Figure 6 . The position of HD 23478 is indicated by the black 
diamond. The evolutionary tracks (black dotted-dashed lines) and 
isochrones (red dotted lines) assume a non-rotating star of solar 
metallicity ||Ekstrom et al.|2012[|. 


scales associated with, for example, stellar rotation. All of 
the uncertainties are reported to be between 0.005 and 0.011 
magnitudes with a mean value of 0.007. Out of the 75 ob¬ 
servations, three report the presence of anomalous measure¬ 
ments between the FAST and NDAG Data Reduction Con¬ 
sortia. Furthermore, one observation reports a magnitude 
that is 2.7{(jHp) larger than (iFp) and is likely an outlier. 
The exclusion of any of the flagged or outlying data points 
had an insignificant effect on the estimated rotational period 
with a maximal deviation of AP < 5 s. 

The period of the Epoch Photometry observations was 
found by assuming a sinusoidal fit and calculating the re¬ 
duced distribution for periods ranging from 0 to 5 days 
with a resolution of 5P ~ 5 sec. This yielded a range of 
minimal-x^ solutions with the global minimum having a 
reduced value of Xred ~ 0.85, occuring at a period of 
Prot = 1.0498 d. The Hipparcos measurements phased at this 
best-fitting period are shown in Fig. along with a subsam¬ 
ple of the full periodogram centered on Prot = 1.0498 d. 

Assuming a 3a confidence limit corresponding to an un¬ 
certainty of 5P = 0.0004 d, we find that only one other pe¬ 
riod has a x^ value below this threshold at P = 1.02454 d 
where the x^ peak is shown in Fig. (bottom) bounded 
by the dashed dark blue lines. The 1.02454 d period does 
not result in a better fit compared with the 1.0498 d pe¬ 
riod; moreover, a similar analysis has previously been ap¬ 
plied to ground-based photometry resulting in a photomet¬ 
ric period of 1.0499 d (for which no uncertainty is reported) 

( Jerzykiewicz|1993 1. We therefore adopt the 1.0498 d period 
along with the ephemeris 


JD = 2448700.606 ± 1.0498(4) • E (1) 

where the reference JD (2448700.606) corresponds to the 
epoch of Hp maximum. 
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Phase 



1.02 1.03 1.04 1.05 1.06 1.07 1.08 

Period [d] 


Figure 7. Top: Hipparcos Epoch Photometry phased with a pe¬ 
riod of 1.0498d. Bottom: A subsample of the full reduced ver¬ 
sus period centered on the period of the best fitting sinusoid. The 
1.0498 d minimal-x^ solution is bounded by the dashed dark blue 
lines. The dashed light blue lines (Prot = 1.02454 d) indicate the 
only other peak falling below our chosen Str threshold. 


When phased with this ephemeris, the Hp data ex¬ 
hibit a roughly sinusoidal variation with a full amplitude 
of 0.015 mag and extrema located at phases of 0.00 ± 0.05 

and 0.50 ± 0.05. _ 

As discussed by Jerzykiewicz (19931, several possibili¬ 
ties exist for the origin of the 1.0498 d photometric period. 
Slowly pulsating B stars typically have periods ~ 0.5 — 5.0 d 
and can be conclusively identified by the presence of mul¬ 
tiperiodicity (Waelkens 19911. Jerzykiewicz (19931 noted 
that g-mode pulsations could potentially explain the ob¬ 
served variability. Another explanation is that HD 23478 is 
an eclipsing binary with a dimmer, less-massive compan¬ 
ion for which no spectral signature is found. [Jerzykiewicz| 
(19931 rejected this possibility arguing that radial velocity 
variations > 80 km s“^ would be exhibited in this scenario. 
The observations analyzed here do not exhibit any variations 
^ 5 km s“^. 

Assuming that the photometric variations are instead 
caused by rotational modulation, the inclination angle of 
the stellar rotation axis, i, can be inferred from the pro¬ 
jected rotational velocity, usini = 140 ± 10 km s”'^, and the 
stellar radius, R = 2.7lQg Rq, resulting in an inclination of 
where the uncertainty is propagated from the 
3a value of 5 Prot = 0.0004 d and estimated errors in R and 
vsini. These values imply a stellar rotational velocity of 
V = 150^“ km s ^ which is consistent with typical values 


for B stars ( 

McNally 19651. Therefore, in agreement with 

Jerzykiewicz 

(19931, we conclude that the observed Hippar- 


cos photometric variations are mostly likely generated by 
the star’s rotation. 


5 MAGNETIC FIELD 

The longitudinal magnetic field, (Pz), is measured through 
Zeeman signatures in the circularly polarized Stokes V spec- 

© 0000 RAS, MNRAS 000, 000-000 



Velocity [km s ’] 


Figure 8 . Examples of He+metal LSD profiles at 3 phases show¬ 
ing Zeeman signatures in Stokes V (top, red), along with associ¬ 
ated Stokes I (bottom, black) and diagnostic null (middle, blue) 
profiles. Vertical dotted lines indicate the integration range. 


trum. In order to maximize the SNR in these observations, 
the Least Squares Deconvolution (LSD) cross-correlation 


procedure (Donati et al. 1997 Kochukhov, Makaganiuk & 


Piskunov 20101 was applied to all ESPaDOnS and Narval 


spectra. The LSD line masks are constructed using spectral 
line lists associated with a specified surface gravity, effec¬ 
tive temperature, detection threshold (i.e. the central depth 
of the weakest lines included in the mask), and microturbu¬ 
lence (umic) obtained from the Vienna Atomic Line Database 















































































































































8 


J. Sikora et al. 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Phase (P = 1.0498(4) d) 


Figure 9. Longitudinal magnetic field measurements phased ac¬ 
cording to Eqn. ^ The curves correspond to the best-fitting sinu¬ 
soids. Top: (Bz) values derived using He-pmetal LSD line profiles. 
Bottom: {Bz) values derived from H/3 Stokes V observations. 
Black circles correspond to measurements obtained by Narval and 
ESPaDOnS while red diamonds are those obtained with dimaPol. 


(VALD2) ( Piskunov et al.|1995 I. In the case of HD 23478, a 
detection threshold of 0.01 and Wmic = 0 were used. Terres¬ 
trial atmospheric absorption appears as deep, narrow telluric 
lines in the Stokes / LSD profiles that can significantly affect 
the value of (Bz); therefore, any lines near spectral regions 
where tellurics were present were removed. 

(Bz) measurements are known to vary to some extent 
depending on which element’s atomic absorption lines are 


used to analyze the Zeeman signatures (e.g. Pyper 1969 


Bohlender et al. 1987 Yakunin et al. 20151. Several LSD 


line masks were therefore created in order to better evaluate 
the strength and variability of the inferred field - one which 
contained only He lines, one in which both H and He lines 
were removed leaving only absorption due to metals, and 
one containing both He and metal lines. 

For each line mask, we adopted normalization values 
of the wavelength, Lande factor, and line depth of 500 nm, 
1.2, and 0.2, respectively. Both the pure He and pure metal 
line masks yielded results similar to the combined He-pmetal 
line mask. This latter mask resulted in the smallest uncer¬ 
tainties. Examples of the Stokes V and I, and of the N 
(circularly polarized intensity, total intensity, and null diag¬ 
nostic) LSD profiles obtained using the He-|-metal line mask 
at 3 phases are shown in Fig. The profiles corresponding 
to phases of 0.140 and 0.817 were obtained by combining 
two consecutive measurements resulting in a moderately in¬ 
creased SNR. 

The three LSD profiles shown in Fig. demonstrate 
that no strong variability is observed in the He and metal 
line Zeeman signatures. We also computed the longitudi¬ 
nal magnetic field from each LSD profile using Eqn. (1) 
of Wade et al. (20001 where an integration range of u G 
(—140,170) km s”’^ was used. The values of (Bz) derived 
from the He-|-metal mask are reported in Table[^ We phased 
these measurements according to the photometric period of 
Eqn. 0- The phased longitudinal field measurements cal¬ 
culated using the He-|-metal ((Bz)He+metai) LSD line mask 
are shown in Fig. 

Measurements of the longitudinal field were also de- 


Table 4. Calculated values of various parameters associated with 
the magnetosphere of HD 23478 as derived from He-|-metal LSD 
profiles and H/3. For those values dependent on B^, which diverges 
as /5 —t 0° and i —t 90°, only lower limits are provided. 



He-^metalsLSD 

H/3 

i [deg] 


oy_io 


10-10 Moyr-l] 


2 q+0.8 


Voo [x 10 ® km s-l] 


1 2 +°-^ 


Rk [R.] 


^•^-1.2 


/3 [deg] 



8±8 

Bd [kG] 

> 12.0 


^ 9.5 

V* 

^ 2.1 X 10 


1.3 X 10® 

Ra [R*] 

> 21.3 


^ 19.0 

Ra/Rk 

^ 4.6 


^ 4.1 


termined using the core of the H/3 line. As with the mea¬ 
surements obtained from LSD profiles, the first moment 
method was applied. We integrated the ESPaDOnS and 
Narval Stokes V and 7 profiles across the full core, from 
—200 to 220 km s“^. Stokes I was normalized to the spec¬ 
tral flux computed from the 10 pixels located immediately 
outside the integration range. The process used to calcu¬ 
late the DAO longitudinal held measurements is outlined in 
Monin et al.| | [^12[ ). The ESPaDOnS and Narval measure- 


ments are listed in Table [T] and the DAO measurements are 
listed in Tabled 

The phased H/3 longitudinal measurements ((H^/h/j) 
show a somewhat larger amplitude variation of 0.6 ± 0.1 kG 
compared with the 0.3±0.1 kG amplitude obtained from the 
He-hmetal LSD profiles (Fig. [^. 

While the sinusoidal fits to the two (Bz) data sets ap¬ 
pear to reproduce the phased measurements reasonably well, 
only the variability observed in {Bz)hi 3 is statistically sig¬ 
nificant. Furthermore, the phases of the extrema of the si¬ 
nusoidal fits (0.4 — 0.5 and 0.9 — 1.0 with uncertainties of 
about 0.05) are in agreement amongst the two curves. This 
suggests that the marginal detections of (Hz/ue-i-metai vari¬ 
ations are in fact real. 

In the context of the Oblique Rotator Model (ORM), 
periodic variations in (Bz) can be understood to be due to 
the rotation of a star having a dipolar magnetic field com¬ 
ponent with an axis of symmetry that is inclined from the 
rotational axis by an angle /3 ( Stibbs|[l950 1. Assuming that 
HD 23478’s magnetic field is characterized by an important 
dipole component, as typically observed in hot stars, the 
faint variability observed in the (Bz) measurements shown 
in Fig. [^implies an obliquity angle of /3 « 0°. Other expla¬ 
nations for the weak longitudinal field variability can be ex¬ 
cluded based on the inferred rotational period and observed 
spectral properties. For instance, regardless of the value of /3, 
no variations in (Bz) would be apparent if i = 0°; however, 
this is inconsistent with the i = 69liJ° calculated in Sec¬ 
tion]^ Weak variability would also be observed if the star’s 
rotational period is on the order of the observing timescale 
of Tabs ^ 30 d. This is inconsistent with HD 23478’s inferred 
1.0498d rotational period and relatively large vsini. 

In the ORM, the obliquity angle is derived from i along 
with the best sinusoidal fit to the (Bz) measurements, given 
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by {Bz) = Bq -\- B\ s.m2Ti[(f) — (j>o), using the relation given 
in Eqn. 3 of Preston (19671. A least-squares analysis yielded 
the fit to {i3z)He-i-metai shown in Fig. described by the 
parameters Bo = —2.1 ± 0.1 kG, Bi = 0.3 ± 0.1 kG, and 
<j)o = —0.3 ± 0.4 where the uncertainties correspond to Icr 
confidence. Applying Eqn. 3 of Preston (19671 with the in- 
clination angle derived above and r = 0.77 ± 0.12 (Eqn. 2 of 


Preston 
value of 


19671 results in an obliquity of /? = 3i'i3°. A similar 
8 ± 8° was found for {Bz)h/ 3 . 

With the inclination and obliquity angles known, the 
polar strength of the surface dipole component of the mag¬ 
netic field can be estimated using the method of [Preston] 
| |1967| ). For this calculation, we use a maximum longitudinal 
field magnitude of = 2.4 ± 0.2 kG as given by the 

sinusoidal fit to (i3z)He+inetai and a limb-darkening coeffi¬ 
cient determined by linearly interpolating the values asso¬ 
ciated with the TLUSTY model atmosphere grid yielding a 
value of u = 0.398 for the Teff = 20 kK, logg = 4.2 model 
i jDaszyhska-Daszkiewicz fc Szewczuk||2011[ ). This yields a 
dipole field strength of Bd ^ 12.0 kG where no upper limit is 
provided due to the divergence of the dipole field as /3 —>■ 0° 
and i —>■ 90°. A relatively small difference in the lower limit 
of Bd was found (9.5 kG) using the combined Narval, ES- 
PaDOnS, and DAO {Bz}hi 3 values. The results of the mag¬ 
netic analysis are summarized in Table 

We note that an upper limit on the dipole field strength 
may be estimated based on the fact that Zeeman splitting is 
not observed in the spectral lines themselves. However, given 
the high usini of the star, those constraints are not likely to 
be very helpful. Moreover, the field strength may be better 
constrained by generating synthetic Stokes V profiles (e.g. 


Petit & Wade 2012 l. On the other hand, its utility would 


likely be hampered by the weak variability of the Stokes V 
profiles. Such an analysis is beyond the scope of this work 
but could be carried out in the future using a more extensive 
and higher SNR spectropolarimetric data set. 


6 EMISSION AND VARIABILITY 


HD 23478 was originally identified as a potential GM star 
based on the large degree of emission occuring at all near- 


infrared H lines (Eikenberry et al. 20141. As expected, the 


most significant emission detected in the visible observa¬ 
tions presented here occurs at the Hq line. Fig. shows 
the minimum and maximum emission observed at phases 
of 0.740 and 0.516 and visible at all phases forming a dou¬ 
ble peak with a separation of approximately ±480 km s“^. 
The large Doppler shifted velocities of the peaks are more 
than 3 times the star’s projected rotational velocity of 
usini = 140 ± 10 km s“^ (indicated by the dotted vertical 


lines in Fig. 101. If we assume that this emission is produced 


by gas bound in corotation with the star, this implies the 
presence of an emitting plasma at R> 377*. 

Although Ha presents the largest emission throughout 
the observed spectra, it varies relatively weakly (in compar¬ 
ison to a Ori E, for example) over the 1.0498 d photometric 
period. This is demonstrated by the equivalent width (EW) 
measurements shown in Fig.[^and the overplotted line pro¬ 
files shown in Fig. |12| Similar to the ORM interpretation of 
a sinusoidally varying longitudinal magnetic held, the con¬ 
tinuous and strong emission of CM stars is also expected 



-1000 -500 0 500 1000 

Velocity [km/s] 


Figure 10. Maximum Ha emission (solid dark blue), as observed 
at a phase of 0.740, forms a double-horned peak separated by a 
velocity of ±480 km s“^ well beyond the stellar surface indicated 
by the vertical dotted black lines at usin* = ±140km s”'^. Min¬ 
imum H« (solid light blue) observed at a phase of 0.516. The 
dashed black curve is the Tgff = 22kK, logg = 4.4 TLUSTY 
model spectrum discussed in Section|3.2| 


to be rotationally modulated (e.g. Oksala et al. 20121. We 
find that the Hq EWs vary coherently when phased by the 
1.0498 d period and is therefore consistent with this inter¬ 
pretation. 

The clearest example of HD 23478’s spectral variability 
is found in He line EW measurements as demonstrated by 
He I A4713 and He i A5876 shown in Fig.[^ While low SNRs 
and sparse phase coverage prevent an independent determi¬ 
nation of Prot from these EW variations, our observations 
are in agreement with a 1.0498 d period. Furthermore, Ha 
and the He lines do not vary coherently when phased by 
the other candidate 1.02454 d rotational period discussed in 
Section This provides further evidence in support of a 
1.0498 d rotational period. 

The maximum and minimum Ha EWs occur at phases 
0.39 ± 0.05 and 0.89 ± 0.05, respectively. These values are 
approximately in phase with both the Hipparcos photome¬ 
try and the (Pz) measurements. The He l A5876 and He l 
A4713 EW measurements are in mutual agreement and 
have maximum and minimum values at 0.76 ± 0.06 and 
0.26 ± 0.06. They are shifted in phase relative to Hp and 
(Pz) by 0.24 ± 0.08, i.e. one-quarter of one cycle. 


7 MAGNETOSPHERE 


The defining characteristic of CM stars aside from the often 
present strong and broad emission lines is their large ratio 


of the Alfven radius to the Kepler co-rotation radius (Petit 


et al. 20131. The Kepler radius defines the point at which 


a rigidly rotating star’s gravitational force is balanced by 
the outward centrifugal force. Based on HD 23478’s inferred 
mass and rotational period, we calculate a Kepler radius of 
Rk = 2.9lj'2 ^0- The Alfven radius {Ra) can be estimated 
from the wind confinement parameter (t;*) which depends 
on the equatorial magnetic field magnitude (Peq), the stel¬ 
lar radius, mass-loss rate in the absence of a magnetic field 
{Mb=q), and terminal wind speed (IGo) as derived by ud 
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Figure 11. Equivalent widths of variable absorption lines phased 
with the star’s rotational period of Prot = 1.0498 d. 
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Figure 12. Overlayed line profiles associated with the EWs 
shown in Fig. from each of the 12 measurements. Significant 
variability was found in all identified He I lines. 


|Doula, Owocki fc Townsend| (|2008[ ) . Following the treatment 
carried out by Petit et al. (20131 . we calculate n* using the 


calibration for Mb=o derived by Vink, de Koter fc Lamers 


( |2000[ ) with Voo = l.SVesc for a Teff = 20kK star. 

Assuming the Bd values calculated in Section]^ we de¬ 
rived ry* and Ra based on the He-l-metal and H/1 analyses 
and find the results to be consistent. The larger uncertain¬ 
ties on the H/3 measurements yielded the lowest minimum 
and Ra at 1.3 x 10® and 19.0 i?*, respectively. With a 
Kepler radius of Rk ~ Rq, we obtain a ratio of 

Ra/Rk = 13.3 (with a lower limit of 4.6) thereby placing 
HD 23478 well within the CM star regime (i.e. those with 
Ra/Rk > !)• The full range of magnetospheric parameters 
calculated from (Bne+inetai) and (Bh/s) are listed in Table 

H 


In the context of all known CM stars, the estimated 
Kepler and Alfven radii place HD 23478 within a region of 
the magnetic confinement diagram - Fig. 3 of [Petit et al.| 


(20131 - populated by approximately 12 mid to early B- 
type stars. Six of these are reported to exhibit Htr emission 
and only one from this subset - HD 142990 - has a rotational 


period shorter than that of HD 23478 (Bychkov, Bychkova 
fc Madej] 2005 (. The lowest reported dipolar magnetic field 


strengths of these 12 CM stars is HD 176582’s Bp Is 7.0 kC, 
a Bp star that exhibits comparatively weaker Ha emission 
( Bohlender fc Monin|2011| . 

We note that HD 23478 exhibits similar physical and 
spectral properties to the well-studied B2Vp star u Ori E. a 
Ori E has a comparable effective temperature of 23 ± 1 kK 
( Groote fc Hunger] 1982 1 and its rotational period and pro¬ 
jected rotational velocity differ from HD 23478 by as little as 
13% and 6%, respectively ( Townsend et al(]|2010 |. Ha EW 
measurements of a Ori E show stronger variability ~ SA 
( Oksala et al.|[2012 l compared with the AEWhq ~ 1 A de¬ 
rived for HD 23478. In terms of the empirical Ra/Rk 10 
limit noted by Shultz et al. (20141 for the occurence of Ha 
emission in CM-hosting stars, these differences in emission 
properties are consistent with a Ori E’s likely higher Rk/Ra 
of approximately 15 (Petit et al. pol^ . 

Compariso ns with the rigidly rotatin g magnetosphere 
(RRM) modeQ ( jTownsend fc Owocki|2005 \ suggest that the 
photometric variability of HD 23478 may be caused by vari¬ 
able occultation of the stellar disc by the magnetosphere. In 
this scenario, the plasma forms a circumstellar disk in the 
magnetic equatorial plane. A non-zero obliquity may then 
allow the plasma to periodically eclipse the stellar disk re¬ 
sulting in an observed dimming of the star. 

The Ha emission and unsigned longitudinal field 
(l(Bz)l) associated with an i = 70° and /3 = 10° RRM 
model are predicted to be in phase; this is consistent with 
our analysis of HD 23478. This supports our interpretation 
of the emission and its variability as due to magnetically con¬ 
fined plasma in a CM. On the other hand, the model also 
predicts a maximum photometric brightness at the phase 
of maximum |{B 2 )|, whereas we observe roughly the oppo¬ 
site. Therefore, the periodic dimming of the star is likely not 
related to magnetospheric occultation. 

Photometric variability may be observed if chemical 
spots are present on the star’s surface, which serve to re¬ 


distribute the star’s flux into the UV (e.g. Peterson 1970 
[Krticka et al.|[^013[ 2015[). The variable He line profiles in 


troduced in Section ^ provide evidence for He spots on the 
surface of the star. However, the He EWs do not vary in 
phase with Hp but are instead shifted by one quarter of 
a cycle. Thus, He spots are also not able to fully account 
for the observed Hp variations. On the other hand, spots of 
other elements (in particular Si and Fe) may well also con¬ 
tribute. Such an investigation would be a useful element of 
a future, more detailed study of HD 23478. 


8 CONCLUSIONS 

Based on our analysis, we draw the following conclusions: 

(i) HD 23478 is a main sequence magnetic He-strong star 
(Teff = 20 ± 2kK) exhibiting peculiar line strengths and 

^ for a visualization, see: 

WWW. astro, wise. edu/_townsend/static.php?ref=rrm-movies 
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distorted line profiles of He, Si, and Fe, indicating likely 
nonuniform surface distributions of these elements; 

(ii) the 1.0498(4) d period found in Hipparcos epoch pho¬ 
tometry is most likely the star’s rotational period which con¬ 


firms the earlier findings of Jerzykiewicz (19931; 


(iii) the presence of strong Ha emission is consistent with 


the discovery by Eikenberry et al. (20141 of emission in Ha 


and nIR H lines. This emission forms a double peaked pro¬ 
file separated by u ~ 960 km s“^ which likely originates from 
plasma located at a distance of i? > 3 R*. We observe sig¬ 
nificant variability in Ha and in various He lines. Measured 
EWs of these lines are coherently phased by Prot = 1.0498 d; 

(iv) the surface magnetic field is inferred to have an im¬ 
portant dipole component, with a minimum polar strength 
9.5 kG, nearly aligned with the star’s rotation axis. 

Our results unambiguously imply the presence of a cen- 
trifugally supported, magnetically confined plasma around 
HD 23478, and therefore confirm the hypothesis of |Eikenf^ 
berry et al. (20141 that it is a member of the growing class 


of magnetic B stars hosting centrifugal magnetospheres. 

Our analysis yielded a relatively significant discrepancy 
in the determination of the star’s effective temperature. lUE 
spectra used in the SED modelling imply Tes = 22 kK while 
the observed Balmer lines require a lower temperature of 
18 kK to provide the best fit. Chemical peculiarities, line 
profile distortions, and relatively weak and scarce metal 
lines prevent this ambiguity from being confidently resolved 
within the scope of the present study. 

The low inferred obliquity angle of /3 ^ 16° paired 
with the high 69l^J° inclination angle introduce basic un¬ 
certainties in modelling the star’s magnetic field and mag¬ 
netosphere. Consequently, only lower limits on the Alfven to 
Kepler radius ratio and the magnetic field’s dipole compo¬ 
nent were obtained. We find that the derived Ra/Rk ^4.1 
and Bd ^ 9.5 kC place HD 23478 within the more extreme 
subset of stars hosting centrifugal magnetospheres. 

Future work is necessary to more precisely constrain 
the magnetic field topology, magnetospheric properties, and 
chemical peculiarities of this star. Although the class of 
strongly Ha emitting stars hosting centrifugal magneto¬ 
spheres is growing, the total number is still relatively low. 
Therefore, a detailed characterization of all known cases is 
important in order to better understand the interaction be¬ 
tween magnetic fields and the winds of hot stars. 
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